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Abstract 



We argue that the RSI model can provide a strong signature in gravitational waves. 
This signal is a relic stochastic background generated during the cosmological phase 
transition from an AdS-Schwarschild phase to the RSI geometry that should occur at 
a temperature in the TeV range. We estimate the amplitude of the signal in terms of 
the parameters of the potential stabilizing the radion and show that over much of the 
parameter region in which the phase transition completes, a signal should be detectable 
at the planned space interferometer, LISA. 



1 Introduction 




The nature of the theory underlying the electroweak phase transition will hopefully be re- 
solved within the next five years at the LHC. However, other indirect probes of the weak 
scale could supplement these results and provide further important insights. The nature 
of the dark matter, for example, might provide weak scale information. In this light, the 
planned LISA gravitational wave detector could be very exciting since the frequencies of the 
observable gravitational waves lie in just the right range for exploring the electroweak scale. 
In this letter, we apply this observation to the warped five- dimensional spacetime of the RSI 
model and show that the early universe phase transition between the AdS-Schwarshild and 
RSI phases could provide a sizable signal. 

The frequency of gravitational waves observed today is 

6 x 10~ 3 mHz ( -^V /6 T * A (i) 
U00/ 100 GeVH* V ; 

where /*, T*, H*, g* are respectively the characteristic frequency, temperature, Hubble fre- 
quency and number of relativistic degrees of freedom at the time when the gravitational 
waves were produced. For weak scale temperatures and ~ 10 2 (as expected for weak 

scale phenomena as we explain below), this is peaked in the LISA band (10~ 4 — 10~ 2 ) Hz 
and is actually a stronger signal for the range slightly outside that probed by the LHC. 
That is, even if weak scale physics is at the high end of the LHC range, we might probe the 
underlying theory further at LISA. 

However, not all weak scale physics is relevant to LISA observations since only rather 
dramatic dynamical phenomena will give rise to detectable gravitational waves. One such 
possibility is a sufficiently strong first order phase transition. In this letter, we point out that 
if the RSI warped geometry resolution to the hierarchy problem [1] is correct, it could yield 
strong detectable signals. This would be significant not only as a potential confirmation of 
LHC results, but is also very promising because the potential energy reach of the gravitational 
wave signal is higher, and phase transitions up to 10 TeV might yield visible signals, while the 
second generation space interferometer Big Bang Observer (BBO) can explore even higher 
energies up to 10 7 GeV [2]. Furthermore, the LISA detector would probe the early universe 
phases, potentially supplementing whatever might be learned at the LHC. 

The key point is that at high temperature, there is an AdS-Schwarschild phase involving 
a single brane where the graviton amplitude is peaked, whereas at low energies, there are 
two branes with a slice of bulk AdS in between [1]. In Ref. [3], it was shown that in the 
perturbative regime one expects a first order phase transition between these two phases, 
which proceeds through the nucleation of "brane bubbles" . Ref. [3] focused on the difficulty 
of completing this phase transition in the perturbative regime consistently with small back- 
reaction, finding their analysis favored lower N (where N is the parameter of the conformal 
field theory that determines the ratio of the five- dimensional Planck scale and the AdS scale). 
In this letter, we discuss some aspects of these results, considering also the negative e case 
(essentially the squared mass of the Goldberger-Wise scalar living in the AdS bulk), and 
demonstrate the potential detectability of the gravitational wave signal associated with the 
first order phase transition over a large region of parameter space. 

Our analysis follows closely the methodology of Refs. [4-11] which applies to very strong 
phase transitions (like the phase transition we are considering in this work, as will be shown 
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shortly), in which case bubble expansion proceeds via detonation [7]. In this regime, the 
gravitational wave signal depends only on two parameters. The first is the dimensionless 
parameter a, which is defined as the ratio of the latent heat to the radiation energy density 
evaluated at the nucleation temperature T n . To achieve a visible signal at LISA, a must be 
sufficiently big, at least 0.2 [2,11]. 

The second parameter, (3, tells the time variation of the bubble nucleation rate and 
hence the length of time of the phase transition. It is defined as (3/H = Td(S%/T) / dT where 
5*3 is the free energy of the bubble and the derivative is also evaluated at the nucleation 
temperature. A large signal requires a relatively slow phase transition so [3/H should be 
small. A visible signal at LISA requires (3/H < 10 3 . (3/H is dimensionless and its size 
is mainly determined by the shape of the effective potential at the nucleation temperature 
T n . It depends on the energy scale T n only logarithmically. Typically, (3/H ~ S3/T. Nu- 
cleation occurs when the probability to nucleate one bubble per horizon time and horizon 
volume ~ T 4 e~ Sz l T / H A becomes of order unity. Therefore the value of S3/T at T n is about 
A\sx(m P i/T n ). For T n ~ 10 2 - 10 3 GeV, we then expect (3/H ~ O(10 2 ) which is in the visible 
range. 

There are two sources of gravitational waves from a first-order phase transition: bubble 
collisions and turbulence in the plasma. The corresponding relic signals, expressed in terms 
of the fraction of the total energy density today, at the peak frequencies f coll and / turb , 
are [7-9,11] : 
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u s and k are respectively the bubble wall velocity, the turbulent fluid velocity and the fraction 
of vacuum energy which goes into kinetic energy of bulk motions of the fluid. They are given 
by [7,11,12] 
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Formulae (j4j), (j3j) for turbulence were recently corrected in [13] and they are being revisited 
for the collision signal [14]. However, for the purpose of the present paper, we use Eqs. (J2]), 

©, dm, ©. 

The phase transition occurs when the radion that determines the distance between the 
two branes of RSI is stabilized. Note that the original description of RSI is five-dimensional, 
yet we are using a four-dimensional formalism to determine the strength of the phase tran- 
sition and the gravitational wave signal. This is justified in the regime where the radion is 
light (lighter than the KK modes) so that it dominates the potential in the RS regime, which 



2 



requires that the parameters determining the radion potential are perturbative. In the high 
temperature phase, the AdS-Schwarschild metric can be interpreted holographically in terms 
of a four- dimensional CFT, so one can use the four-dimensional formalism in that regime as 
well. 

To evaluate the relevant quantities, we use the formalism of Creminelli et al [3], who 
assumed the original Goldberger-Wise (GW) model in which there is a scalar field whose 
bulk mass m is determined by the parameter e = ^4 + m 2 /k 2 — 2. This field takes values vo 
on the Planck brane and V\ on the TeV brane. From the five-dimensional perspective, the 
radion is stabilized by the tension between the mass term and the gradient contribution to 
the energy. Because of the background anti de Sitter space, e ~ m 2 /4k 2 can be positive or 
negative, though the original paper considered only the positive case [15]. Other variations 
on this model include interactions of the scalar field but we don't expect this to change the 
results significantly. 

One can also consider the holographic interpretation of this model [16,17]. In this in- 
terpretation one has a marginal operator, which is relevant for positive e and irrelevant for 
negative e. The positive case is peculiar from a holographic perspective. In this interpreta- 
tion the breaking of conformal symmetry in the IR occurs because of the competing effect 
of two terms. The negative e case is more general and conventional, and corresponds to a 
coupling that grows in the IR, where it breaks the conformal symmetry, analogously to the 
QCD phase transition. We will consider both positive and negative e. 

The five- dimensional metric is ds 2 = e~ 2ky ri /iu dx tl dx 1/ + dy 2 where k = 1/L is the AdS 
curvature. The kinetic term for the radion field /i = ke~ ky , in terms of the five- dimensional 
Planck mass M, is [17-19] 

C kin = -12v^(ML) 3 (9/i) 2 (6) 
and its induced four-dimensional potential is [15] 



V GVf ( f i) = ev^ + (4 + 2e)/i>!- V/V/^T-o^ + ^/i 4 +0(jj»/f4) 



(7) 



where /xq is the UV scale and |e| <C 1 has been assumed. The terms Vq and V\ are the vevs 
in Planck units of the Goldberger-Wise field on the Planck and TeV branes. 

We sometimes modify the potential of [15] to allow a term <5Ti/i 4 corresponding to a 
change of the TeV brane tension. Such a term is permissible and allows for a larger range 
of viable models. 

Provided that STi < evf the potential above has a global minimum at 



/^TeV± 
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where ± corresponds to the cases e > and e < respectively. For e < 0, we have 
the additional condition STi > —vf(A + e) so that the minimum above does not become a 
maximum. For e > and 8Ti = 0, /i To v± ~ Ho i v i/ v o) 1 ^ ■ 111 that case, the hierarchy between 
the weak and Planck scale can be naturally obtained for parameters not far from 1. For e < 
and 8T\ = the only minimum of the potential is at \x = 0. With nonzero 8T\ < there is 
also the desired minimum at fi ~ ^o( v o/ Vi) 1 ^^^ — y/8Ti/2vi) 1 ^^ and a viable solution can 
be obtained for vq < Vi [16, 17]. Note that /j = is also a minimum of the potential but for 
e > and 8Ti < — uf(4 + e), the barrier disappears and \i = becomes a maximum. 
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For small e, the value of the potential at the minimum (jSJ) where the radion achieves its 
vacuum expectation value is 

V± « /4 V± e((JTi/2 =F v 1 ^-8T 1 + ev 2 1 ) (9) 

where we have dropped ei^/i 4 , which is common to the high and low energy phases and won't 
contribute to a or j3. 

When 8T\ = 0, e has to be positive. The value of the potential at the minimum scales 
like — e 3 / 2 vj[i* eV , where V\ is the assumed value of the Goldberger-Wise field on the TeV 
brane. As the free energy of a critical bubble is smaller if the minimum is deep, it will 
clearly be smaller when v\ is larger. However, overly large v± could result in a large back- 
reaction to the potential as we will see in Section SJ A quick way to extrapolate the 8T\ = 
results for large 8T\ is to note that the potential at the minimum is then proportional to 
/j,1 CeV e(5T 1 /2^fv 1 y/—5Ti). In other words, the quantity e 3 ^ 2 vj which appeared in the minimum 
energy for 8T\ = is replaced by e8T]v\. Clearly, larger \8T\\ yields a deeper minimum. This 
makes the tunneling amplitude bigger. 



2 Completion of the phase transition 

The free energy of the AdS-S phase is [3] 

F AdS _ s = -2vr 4 (ML) 3 T 4 (10) 



(flOl) is a minimum of the free energy of AdS-S corresponding to the Hawking temperature 
at the horizon of the black hole solution equal to the temperature of the universe. By 
holography, F AdS _ s can be interpreted as the free energy of a strongly coupled large N CFT 
with N 2 = 167r 2 (ML) 3 . The exact relation depends on the precise theory; this formula 
is for M = 4 SU(N) super Yang Mills [20]. We often choose to use the parameter N to 
characterize the AdS curvature. 

In the Randall-Sundrum phase at high temperature, the TeV brane is pushed by thermal 
effects to the AdS horizon in the absence of the Goldberger-Wise field. Since it costs energy 
for n to be nonzero, there exists an energy barrier which leads to a first-order phase transition. 
The phase transition can take place after the energy of the AdS-S phase equals that of the 
RS phase. Another argument for the first order phase transition is provided by the AdS-CFT 
correspondence that relates the RS model with a 4D confining gauge theory. It is well-known 
that the confining phase transition of large N (N > 3) gauge theories is first order (growing 
more strongly as N is large). In the 4D description, the stabilized radion is some glueball 
state. The critical temperature is the temperature at which the RS energy at the minimum 
of the GW potential equals the AdS-Schwarschild solution energy. We have 

T c= (11) 



7T 



'■N 2 J 



The phase transition can proceed only if the bubble nucleation rate is larger than the expan- 
sion rate of the universe i.e. S$/T < Aln(mpi/T n ) where £3 is the free energy of a critical 
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bubble. For electroweak scale temperatures, this corresponds to the condition S 3 /T < 140. 
The latent heat is of the order of V±. 

We now follow the assumption of [3] , where the contribution from the AdS-S side to the 
thermal bounce action is neglected. This is reasonable when v\ and e are small and the 
potential is consequently shallow, in which case the bubbles are big and most of the action 
comes from the RS regime as \x changes from to fiTeV- We therefore integrate only over 
the RS side when evaluating the action. This approximation is reasonable provided that 
T c <C /i To v± i n which case only the radion mode contributes to the action and is valid so long 
as V\ is smalfl- 

The quantity a depends on the free energy difference between the two minima, which is 

V± \T C , 

where we have deliberately left open the possibility that at nucleation Vaw(^) is n °t exactly 
the value of the potential at the minimum V±. 

The quantity j3 is determined from the bubble action, which is greatly simplified in two 
limiting cases: the thin-wall and thick-wall approximations. In the first case, the bubble 
radius is much larger than the thickness of the wall (region over which the value of \i varies). 
This applies when AV± is much smaller than the height of the barrier separating the two 
minima (which we do not know). 

However, if the thin-wall action is too big to allow nucleation, the action proceeds via 
thick-wall bubbles. As the temperature drops, AV± increases and eventually becomes larger 
than the barrier. It is then favorable to make the wall thickness comparable to the bubble 
size to minimize the surface term. We will soon see that the phase transition occurs primarily 
in the thick- wall regime in which supercooling is relatively large. 

The tunneling rate is proportional to exp(— S), where S = S3/T if thermal bubbles 
dominate, where S3 = An J r 2 dr[(dfi/dr) 2 /2 + AV(fi)} is the 3D Euclidean action of an 
0(3)-symmetric critical bubble, and S = S 4 = 2n 2 J r 3 dr[(dfi/dr) 2 /2 + AV(/i)] in the regime 
where 0(4)-symmetric bubbles dominate. 

In general, nucleation proceeds via thermal bubbles but if the phase transition doesn't 
complete before a temperature T such that T < (2i?) _1 , where R is the radius of the 0(4) 
bubble, and S4 < S3/T the 0(4) symmetric solution [22] applies. We find that S4 is never 
smaller than S3/T in the thin wall approximation. However, for thick wall, there will be 
some region of parameter space where S4 is more favorable. 

In the thin-wall approximation [22], 

S 3 = (16rr/3) (3iV 2 /27r 2 ) 3/2 S 3 /(AV ± (/i)) 2 (13) 

where Si = Jo TeV± d^2\ AV\ ~ — fo TcV± dfi^—2V± is the surface tension evaluated in the 
limit T -> T c and the (3A^ 2 /2tt 2 ) 3 / 2 factor comes from the canonical normalization of fi. We 
evaluate the denominator in (fl~3l) at the minimum of the potential, fi = fi TcV , for which: 

AT4 = 14(1 - (T/T c ) 4 ) (14) 
This leads tcfl S 3 ~ (167r/3)2 3 / 2 (/i 3 oV /w^g) (3A^ 2 /2vr 2 ) 3/2 /(l - (T/T c ) 4 ) 2 which can be ex- 



1 Before the holographic description was emphasized, a study of the phase transition with a different 
radion potential (and a different initial thermal phase) was performed in [21]. 
2 This is a factor 8 x 2 3 / 2 larger than in [3]. 
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pressed as 

S 3 ^ 2.95 N 7 / 2 (TJT) (i5) 



T 



3/4 (1 _ (T/Tc 



4\2 



In this formula, we have neglected the unknown barrier contribution which could change the 
result by a factor of order unity. It turns out the transition is dominated by the thick wall 
regime even without additional suppression so that this formula is adequate. 

The T-dependent factor in the right of eq. (j!5p blows up as T approaches T c . As T 
slightly decreases, it becomes of order one, reaches a minimum at T = Te/3 1 / 2 and grows 
again. Clearly, the nucleation condition S 3 /T < 140 is satisfied when this factor is near its 
minimum, for which it is ~ 0(1). This typically takes place when T is between Tc/3 1 / 2 and 
0.9 T c . 

From eq. f)15p . it is clear that there is a constraint that N not be too large (or equivalently, 
the AdS curvature scale not be too small). This follows from the large entropy in the AdS- 
Schwarschild phase, which grows as N 2 , whereas the RS phase entropy does not. This means 
that the phase transition is entropically disfavored for large N. We can understand the TV- 
dependence in the action as arising from the radion-normalization and from the shallowness 
of the potential. This latter ^-dependence comes from the size of the bubbles, which grows 
with N, and would not be present if V\ was as large as N. However perturbativity of the 
model requires v\ < N. We will comment further on this point below. 

We have ignored the T-dependent corrections to the radion potential so far. We can 
model the radion potential at finite temperature just assuming that the mass of the CFT 
degrees of freedom is proportional to the vev of fi, m — gfi, according to the general formula: 

at// ™n ^ 9fT 4 ^ (-1)" m F 2 nm F g B T 4 ^ 1 m B 2 nm B 

AV(/i ' T) = ^ "2^ ? K2( — } - ? ~W ? n 2 ^ K2i — ] (16) 

This is the 1-loop thermal corrections to the radion potential taking into account the inter- 
actions of the radion with the iV 2 CFT degrees of freedom. J2b/f9b/f = 45iV 2 /87r 2 . At 
/i = 0, this reproduces the T 4 radiation energy density, but as soon as /i is of order T, the 
CFT degrees of freedom are massive and the T-dependent piece falls off. As illustrated on 
the fourth plot of Fig. [3l at T close to T c , the barrier is big and prevents tunneling. But 
as T goes below T c /2, the barrier becomes much smaller than the energy difference between 
the two minima and we tend toward the zero temperature radion potential. The transition 
does not take place until the temperature is low enough so that we are in the thick wall 
approximation. This trend is general, even though detailed predictions will depend on the 
precise coefficient g between the mass of the CFT degrees of freedom and the radion. We 
took g — 1 in the right bottom plot of Fig. [31 

Thin-wall bubbles generally yield too big an action for nucleation and we now show that 
thick-wall bubbles dominate. As we have just argued, it is a good approximation to use the 
zero-temperature radion potential to estimate the tunneling action. The action for thick- 
walled bubbles is S 3 ~ 2-wRy 2 - 4tt R 3 \AV \/3 [23] and its minimal value S 3 = SnR 3 \AV\/3 
associated with R = n^3N 2 /2% 2 / yj2\ AV\ leads to 

\ 3/2 



S 3 « (47r/3)(/i 3 /y / 2]AKj) (3iV 2 /27r 2 ) (17) 
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which is, assuming again that at nucleation \x = /i To v : 

S 3 « (47r/3)( /UTcV 3 / v /2T^) (3iV 2 /2vr 2 ) 3/2 /(l - (T/T c ) 4 ) 2 (18) 

This action is smaller by an overall numerical factor (sixteen) than the thin wall formula. 

We also consider the possibility of 0(4)-symmetric vacuum bubbles. These can be rel- 
evant only if the nucleation does not take place via thermal bubbles. The thin-wall action 
in this case is [22] 5 4 = 27n 2 Sf x (3A^ 2 /(2tt 2 )) 2 /2| AV±| 3 whereas the thick-wall action is 
Si = 7r 2 i?V - tt 2 R a \AV\/2, which is minimal for R = ^3N 2 /2ir 2 / ^\AV\ leading to 

S A = 7r 2 fi\3N 2 /27r 2 ) 2 /(2\AV ± \). (19) 

and if /i = fi TcV : 

S, = n 2 ^ TcV (3N 2 /2n 2 ) 2 /(2\V ± \(l - a 4 )). (20) 

The ratio of the thick- wall to thin- wall action is ^-independent and for \x = fi TcV , it scales 
as: 

S, _ 3V3 I (2 + e) 1 / 4 

S 3 /T yjl - (J) 4 Vbe^[-c T ij4c + e(4 + e)]W 

It will thus be more favorable to tunnel via 0(4)-symmetric bubbles if T/T c is small, e is large 
and b = Vi/N and c = |<5Ti|/u 2 are large. Using the 0(4) rather than the 0(3) symmetric 
solution corresponds to T n < (2i?) _1 . 

The value of /i to be used in eq. ffT7|) and f|T9l is actually not fi TcV but the value at 
the time of the tunneling which is typically smaller than /i Tc v That is because the value 
of the field that minimizes the action is not exactly the value of the field at the potential 
minimum [23]. The field tunnels to a value near the minimum and after tunneling rolls to 
the minimum. We first evaluated the nucleation action assuming that \x = /i TcV at the time 
of nucleation, to get some analytical dependence on parameters. To derive our results, we 
do not apply formula (TT8T) and (1201) where the radion value is taken at the minimum of the 
potential. We search for the value of the radion field fi minimizing the bubble free energy, 
which is, for the 0(3) action, 

I AW 



and for the 0(4) action 



where 



and 



" = W (22) 



|AV ± | = 14[a"-^] , =£ (24) 



V± = /4v K 2^ [c ± V 4c + <4 + 6)] (25) 

is the value of the potential at the minimum of the potential and we remind that c = \5Ti/v\\. 
We re-express everything in terms of the variable 

Y={—) e (26) 
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Defining 

1 + | ± Wc + e(l + f) 
X ± = 7 1L (27) 

1 ~r 2 

we get 

V ± (Y) = 4 eV yt t, 2 [(4 + 2e)(l - X ± F) 2 - e - c] (28) 
Eq (l22i) and (123|) can be rewritten as 

2(2 - e)X 2 ± Y 2 - 8X±Y + 4 - e - c - 3a 4 y^e(2 - e) [-c T \/ 4c + e ( 4 + e )] = ( 29 ) 



(2 + e) (4 + 2e)X±y (1 - X±Y) - 2a 4 Y~~ [-c =p ^/4c + e(4 + e)] = (30) 
We then look for the regions in (e, vi) plane where the phase transition completes. We 
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Figure 1: Lines delimiting the region in (e, V\) parameter space where the phase transition 
completes (we have to be above the dashed line for nucleation to take place). The upper line 
comes from using the approximate eq (I2UI) while the lower line results from solving eq. 



evaluate S 3 /T and S4 at the Y solution of eq (129]) and (1301 and follow their evolution 
with temperature. The nucleation temperature is defined by the condition S3/T = 140 or 
S4 = 140. Whether S3/T or S4 first reaches this critical value determines whether 0(3) 
or 0(4) bubbles are nucleated. Note that the 0(4) nucleation region is simply given by 

V\/N > 3N/ (4^^70(0 + e)) where the right-hand side corresponds to the limit T n — > 0. To 
illustrate how much difference this procedure makes, we compare the approximate solution 
(using /i = /i ToV ) and the accurate one in Fig. [T]for some particular choice of parameters in 
the 0(4) case. We get the same kind of discrepancy for 0(3). The biggest effect is for small 
e in which case the potential at the minimum is very shallow. Consequently, it does not cost 
much potential energy to tunnel to the wrong place, namely to \i < fi TeV . 

We have checked in various regions of parameter space that the thick wall approximation 
gives an answer that is indeed close to the exact bounce computed numerically. This is 
illustrated in Fig. [2] and Fig. [3l 

Finally, we note that there exists the possibility of an inflationary phase before the 
phase transition occurs. This is due to the vacuum energy — V± in the AdS-Schwarschild 
phase if we set the cosmological constant to in the RS phase. Inflation could in principle 
start when the temperature is of order T c . However, even if inflation were to occur, the 
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Figure 2: Comparison of the thin and thick wall approximations (dotted lines) with the 
exact solutions obtained by solving for the bounce numerically (solid lines). 




Figure 3: Top: comparison between thick wall, thin wall and exact solutions at fixed e and 
5Ti\ bottom left: exact results for different values of N. bottom right: Typical evolution of 
the radion potential with temperature. The height of the barrier falls off as T goes down. 
For T below T c /2, it is a very good approximation to use the zero temperature potential to 
compute the bounce. 
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inflationary epoch would last only a very short time as the number of e-foldings would be 
only ln(T c /T n ). Therefore, even if the universe inflates above the nucleation temperature, 
the transition would typically complete in less than a single e-folding. This is true even if 
the universe settles into the wrong (/i = 0) minimum at early times. In fact, because we 
neglected the AdS-S contribution to the bounce, the action associated with tunneling from 
this false minimum would be precisely the same as what we've calculated. 

To summarize this section, we have extended the Creminelli et al calculation to the e < 
and STi ^ cases. In addition, we considered the nucleation of thick-walled bubbles rather 
than thin-walled ones, which applies in the regime of large supercooling. We computed the 
free energy of a critical bubble by searching for the value of the radion field that minimizes 
the action, rather than assuming that the radion value sits at the minimum of the potential 
at the time of nucleation. Also relevant in the regime of large supercooling is the nucleation 
of 0(4) symmetric bubbles. These effects improve the nucleation probability. 



3 Gravitational wave signal 

We can now determine a and (3 and hence the gravitational wave signal. We repeat that 
we are working under the assumption that the phase transition is very strong (as will be 
justified below), which is why predictions can be given as functions of a and (3/H* only. The 
signal obviously grows with a as the more latent heat that is released, the more measurable 
is the phase transition. The quantity a can be estimated as 

\AV±\ 1 V±(Y) 
vr 2 iV 2 T 4 /8 a 4 V ± [ ' 

To get a feeling for the results, we will first assume that at nucleation /x = fi TcV , in which 
case AV± = — 7r 2 iV 2 T c 4 (l — (T n /T c ) 4 )/8. Hence we can determine a simply in terms of the 
ratio of the nucleation to the critical temperature as 

* = {Tt-n)m (32) 

The size of the signal also increases with the duration of the phase transition (3" 1 given in 
terms of a = T n /T c by: 

P_ f S 3 (T n )/T n x ^ « 140 x ^ for 0(3) solution 

H \ S 4 (T n ) x ^ « 140 x for 0(4) solution 1 ' 



Interestingly, a and (3/H do not depend explicitly on the parameters of the GW potential, 
e, V\, and STi but only on a = T n /T c . Of course, T n and T c implicitly depend on e, v\, and 
5Ti. Again, we have not included the extra contribution from the T-dependence in the RS 
potential. This would lower the value of the latent heat but also lower the value of T c and 
we expect that the effect on T n /T c and thus a should not be significant. It could change the 
calculation of /3/H however. 

From Figure HJ it is clear that a can be larger than 1 in much of the parameter regime. 
In particular, it gets very large if T n /T c < 1/3 1//4 , when 0(4) symmetric bubbles have to be 
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a yS/H 




Figure 4: a is the ratio of the latent heat to the radiation energy density in the CFT phase 
at the time of nucleation, given by Eq. (132]) . It increases as the ratio of the nucleation 
temperature T n to the critical temperature T c decreases. Second plot is j3/H from Eq. fl33l) 
where can be understood as the duration of the phase transition. The amplitude of the 
gravitational wave signal increases with a and decreases as (f3/H)~ 2 . 



considered. This justifies the approximation that bubble expansion proceeds via detonation. 
We expect large signals at LISA, whose sensitivity requires at least a > 0.2 and (3 /H < 1000 
for observability of the signal [2, 11] (for a general detectability analysis at LISA and BBO 
in the (ct,/3/H) plane, for any given temperature, see [2]). 

All predictions are functions of the nucleation temperature T n which can be computed 
in the e, v± plane. We will present in Fig. [7] the region of this plane in which the phase 
transition can take place. However, there will be strong constraints from perturbativity and 
back-reaction which we discuss in the following section, and that will reduce the region of 
parameter space where predictions of large signals are reliable. We address these constraints 
in the next section. 

We now show some typical values of the quantities that are relevant to the gravitational 
wave signal as a function of V\. Fig. [5] shows the values of T n /T c , T n //i ToV , a, (3/H as well 
as the characteristic quantities of the gravitational wave signal obtained for a benchmark 
point e = —0.25, \8T\\ = vf/2. The predicted gravitational wave spectrum for two particular 
values of V\ are presented in Fig. [HI In this region of parameter space there is a huge signal of 
gravitational waves. While relatively large values oivi/N are needed for the phase transition 
to take place, once we are in the region where there is a phase transition, the smallest allowed 
values of vi/N lead to the largest signals (as much as three orders of magnitude above LISA's 
sensitivity) as they correspond to the largest amount of supercooling. The largest values of V\ 
lead to larger characteristic frequencies no more observable at LISA but within the sensitivity 
range of BBO. These larger values of V\ are worse for perturbativity in any case. Note also 
that in this case, the temperature of the transition is not far from the KK scale given by /i Tc y 
For the particular region of parameters of Fig. [51 nucleation takes place via 0(4) bubbles for 
Vi/N < 1.1 and via 0(3) bubbles for larger V\/N. Since S$/T and S4 scale as a positive 
power of N, large N suppresses nucleation and we present results for the maximum value of 
N, N = 12, beyond which the phase transition cannot complete. Relatively low values of N 
are reasonable to explore since if N were much larger, one risks losing asymptotic freedom. 
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Figure 5: a and (3/H determine the spectrum of gravitational waves. They are calculated 
for some benchmark point: e = —0.25, N = 12, 8T\ = — 0.5vf. The smallest values of V\/N 
correspond to a large amount of supercooling i.e. a small value of the ratio T n /T c . This ratio 
varies between 0.23 for V\/N = 0.7 to 0.87 for v\/N = 1.5. Assuming that /iTcv = 

5 TeV, 

the corresponding nucleation temperatures are in the range 490 GeV - 2700 GeV. We also 
show the peak frequency / peak of the gravitational wave signal from turbulence and f2 peak /t 2 . 
The peak frequency depends on /i TcV , while a, (3/H and f2 pcak /i 2 do not. As shown in Fig. [HI 
this can lead to a spectacular signal at LISA and/or BBO. 



In any case, the size of the gravitational wave signal does not depend explicitly on N but 
only on the amount of supercooling, which means that we prefer to live close to the limits 
of the region in the (e, V\/N) parameter space where the transition completes. While N is 
very important to determine the region of parameter space where the phase transition takes 
place, once N is fixed, the size of the GW signal depends only on how far we are from the 
limits of the underlying region, which depend solely on V\/N and e. And different values of 
N can lead to the same size for the signal. Therefore, the choice N = 12 in Fig. [5] and Fig. [6] 



12 




Figure 6: Relic energy density in gravitational waves produced at the phase transition, 
Qawh 2 , versus frequency for two values of V\ corresponding to the benchmark point of 
Fig. [5j The upper dashed orange line is the LISA sensitivity. The lower one is the sensitivity 
for the second generation Big Bang Observatory. For each set of spectra, the first peak is 
from turbulence while the second peak is from bubble collisions. v\ = 0.7N leads to a very 
strong phase transition (a = 26, j3/H = 21) at a temperature T n = 490 GeV and a huge 
signal at LISA. For larger v±, (3 / H is higher and the peak frequencies are shifted to larger 
values. At some point, the signal is no more visible by LISA but still observable at BBO. 
The v\ = 1.17V spectrum corresponds to a = 1.6, (3/H = 350 and a nucleation temperature 
T n = 1830 GeV. The horizontal lines are gravitational wave spectra expected from inflation, 
for comparison (BBO is mainly planned to detect these gravitational waves), for two different 
scales of inflation. The dashed red curve is the expected irreducible background due to white 
dwarf binaries [26] . 



does not really matter. 

Note that the results depend on the scale /i ToV only through the frequency of the signal. 
Indeed, a does not depend on fi TeV . S 3 /T, and (3/H do not depend on /i TcV explicitly. It 
is only when they are evaluated at the nucleation temperature that some weak logarithmic 
dependence on the energy scale of the phase transition appears. Therefore, the amplitude 
of the signal Vth 2 virtually does not depend on fi TeV . On the other hand, the temperature 
of the transition and thus the peak frequency of the signal is proportional to /i TeV . In the 
particular example of Fig. [51 we plot the peak frequency in the three 3,5,10 
TeV. This shows that it is possible to see the peak of the turbulence signal at LISA for ji TaY 
as large as ~10 TeV. This does not mean that higher values of fi TaV are not accessible. In 
fact, for some e, V\ values, it could still be possible to see at LISA the low frequency tail of 
the signal for /i TcV as large as a few hundreds of TeV. This can easily be seen by translating 
the blue peak of the first plot of Fig. [6] to the right by a factor /i XoV /5 TeV. 
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4 Perturbativity Constraints 



We have seen in the last section that the first order phase transition from AdS-Schwarschild 
to RSI geometry could lead to spectacular gravitational wave signals if the first order phase 
transition completes at about the TeV scale. In this section, we consider the constraints 
from perturbativity and the limits on our perturbative analysis. We express our answers in 
terms of N of a CFT but that is directly related to the ratio of mass scales through our 
assumed relation (ML) 3 = N 2 /16tt 2 . In order to trust the background AdS solution, we 
impose 1/(1Qtt 2 (ML) 3 ) < 1, which leads to the weak bound N > 1. Most phenomenological 
models have relatively small values of N. As we will see, for N > 12, the transition cannot 
complete. Small Vq and small e, which are needed for perturbativity, suppress the tunneling 
rate. We now consider the constraints on these parameters. 

Essentially there is a trade-off between good ranges of each of the parameters. Smaller 
N leads to better values of Vq and V\. The constraints on Vq and V\ come from avoiding too 
large a back-reaction to the AdS energy of the five- dimensional theory, as well as imposing 
the radion- dominance assumption that was critical to the four-dimensional analysis. When a 
vev is too big, there can be large corrections to the GW potential from KK modes and large 
corrections to the radion kinetic term due to mixing with the GW scalar. According to [24], 
the last two sources can be acceptable however. That is, when the full mass matrix and 
kinetic matrix are diagonalized, one can still find dominance of the light radion. Of course, 
in that case, we are not guaranteed that the radion potential takes the form we assumed 
and we are not guaranteed the transition remains first order. However, we see no reason to 
assume that in all cases the first order phase transition would disappear, but of course we 
don't know for sure without a more complete analysis in the large vev regime. 

Here we will focus on the constraint on small back-reaction to the energy. Notice that 
even this is potentially stronger than is necessary. Imposing the constraint Vq,V\ < iV/(47r) 
i.e. requiring that the vevs of the Goldberger-Wise field are smaller than M, the 5d Planck 
scale, clearly suppresses the back-reaction. As [25] have shown, small back-reaction to the 
energy is not necessarily essential to maintaining the hierarchjo Nonetheless, it would be 
best to have small back-reaction so that the original AdS space analysis can be trusted. We 
view these plots as indicative of the proximity of the perturbativity limits. 

The stress tensor for the Goldberger-Wise field <p is 

TfiN = -9M N [-(d<f>) 2 /2 - m 2 2 /2] - d M <f>d N <t> (34) 
Using the metric ds 2 = (dx^dx^ + dz 2 )/(k 2 z 2 ) and <fi(z) = Az A+e + Bz~ e where 

A ~ Zl k i-(^M) 4+2e ' i-(z M) 4+2e 1 ] 

We find 

T%{z) = 2 V/ L4V +2e (4 + 36) + eB 2 z- 2 - 2 *} (36) 

When we compare this with the energy momentum tensor due to the bulk cosmological 
constant A 5 = -2AM 3 k 2 

= (kz)-\ u A 5 (37) 



3 Note that their model with large back-reaction would also have a large vev that could induce corrections 
to the potential that were not discussed. 
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we find the condition 



v 2 < 12(ML) 3 /|e| i.e. ^< 




) 



c 



(38) 



on the Planck brane and 




|(4 + 3e)(l-X) 2 + eX 2 



12(ML) 3 



i.e. 




y / |(4 + 3e)(l-X) 2 + eX 2 



1 



(39) 




on the TeV brane, where X = 1 + | ± \\j—^r + e + e IP> ^] appears in the relation 

between v and ui, Eq. ([B]). 

If instead of evaluating the constraint on the TeV and Planck branes we integrate it over 
z, we get the same as eq. (138]) . For \5Ti\ = 0, X = 1 and the conditions fl38|) and fl39l) 
respectively become w 2 < 12(ML) 3 /|e| and ^ < 12(ML) 3 /|e| but for |5Ti| = 4v 2 , the back 
reaction on the IR brane leads to the stronger constraint 



Notice these constraints are readily understood when 8T\ is zero. If v \ or Vq (normalized 
to a dimensionless quantity through factors of k) are too large compared to the Planck 
scale, the expansion in powers of the Goldberger-Wise field would be invalid. Furthermore, 
if the Goldberger-Wise field enters only through the kinetic and mass terms, we have the 
approximate condition m 2 2 /2 ~ 2k 2 e<f) 2 < 24k 2 M 3 , we get the conditions above for zero 
8T\. With the 8T\ turned on and large, it already implies some back-reaction near the TeV 
brane. 

We plot these constraints in Fig. [7J For positive e, the IR back reaction constraint is 
too big and nucleation never takes place in a regime where we can ignore the back reaction. 
Because positive e requires vq > v\, the UV backreaction condition is usually stronger. 
However, from a strictly phenomenological perspective, this constraint is not essential as it 
comes from the requirement that the RS model applies all the way to the Planck scale. We 
can relax the constraint with a lower cutoff and assume some unknown ultraviolet physics 
beyond that scale. In Figure [7J , we show how the viable region of parameter space in the 
e, V\ plane increases as we lower the UV scale. However, even a UV scale at 10 5 GeV is not 
enough to make nucleation in the perturbative region. 

4.1 Negative e 

We find that results are more favorable for e < 0. This is not because nucleation is easier 
but because the back reaction constraints are weaker. The point is that the constraints on 
v and vi are parametrically comparable. But for e < 0, we have v < V\ (see Eq. [8]), so 
the perturbativity constraint on vq is generally satisfied once the constraint on v\ is metQ. 

4 A bound on vq comparable to (|38[) is obtained in the particular example of Ref. [24], where the back- 
reaction parameter is defined as I — Attvq/N. To guarantee that the radion and KK modes are heavy 
enough and not too strongly coupled to SM fields on the TeV brane, they need I < 10 i.e. v /N < 1. 
Since V\ ~ fo^Tev/Mo) 6 ! this constraint is extremely strong for e > and incompatible with our nucleation 
condition but for e < 0, it is readily satisfied. 



vl < 3(ML) 3 



(40) 
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Figure 7: e > (top) and e < (bottom). The dashed black line delimits the region 
where nucleation of 0(4) symmetric thick- walled bubbles is possible. We have to be above 
this line for 0(4) nucleation to take place. The nucleation contour for thick- walled 0(3)- 
symmetric bubbles is typically slightly above (see the black solid line on the third plot). 
Close to the nucleation line, the phase transition is very strong. The pink lines delimit the 
region where back- reaction is important. Regions that satisfy the back-reaction constraint 
are below these lines. Almost no region remains for e > 0, but with e < 0, some regions 
satisfy the constraints. We have shown how decreasing N enlarges the parameter space. 



As for the IR constraint (1391) . we find regions of parameter space where it is relaxed, for 
instance if STi = 2evf where the right hand side of the inequality (1391) blows up to infinity. 
This follows from the right hand side of eq (1381) which blows up at large |e| as (/io/^Tov)' 6 '- 
All this is illustrated in Fig. [7J Of course we are ignoring higher order terms so we take this 
result as simply indicative of the fact that although nucleation takes place at or above the 
perturbativity limit in general for N > 12, the leading term and higher order terms might 
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conspire to be small. 

Notice that the position of the (pink) lines delimiting the perturbativity regions are N- 
independent while the nucleation (black) lines get shifted to lower values of v\/N at smaller 
N. For instance, for N ~ 3,4, it is possible to obtain nucleation at Vi/N ~ 0.1 rather 
than at v±/N ~ 1 and therefore to enlarge the region of parameter space where the phase 
transition completes in the regime of small back reaction, as illustrated on Fig. [7J 

There is, however, an additional constraint for the negative e the operator that 

breaks conformal invariance gets strong in the IR. Below the scale of strong coupling A, the 
conformal picture is spoiled so perturbativity also requires that the nucleation temperature 
exceeds A, where A should correspond to the value of \x at which v\ takes its maximum value 
consistent with perturbativity. We find this constraint essentially coincides with Vi/N ~ 1. 
Since the strong coupling scale is not a precisely defined quantity, we would hope that 
strong coupling effects are not big before nucleation. The first order phase transition might 
be possible, even with strong coupling effects included. However, if strong coupling effects 
are important, we can no longer be confident about the order of the phase transition. 

5 Conclusion 

There can be a significant signal for gravitational waves from the phase transition from 
AdS-Schwarschild to RSI. Predictions vary by orders of magnitude depending on the region 
of parameter space associated with the scalar potential stabilizing the radion but if the 
transition completes, the signal is likely to be significant. The KK scale that can be probed 
is much higher than at colliders. An IR scale as large as a few tens of TeV (even hundred) is in 
principle reachable at LISA, depending on the region of parameter space. The uncertainties 
in the computation come from the unknown temperature-dependence of the potential and 
the unknown back-reaction effects. 

We have focused on the gravitational wave signal that is indicative of any first order phase 
transition; we have not included any features peculiar to this theory. That means that even if 
this signal is measured, we cannot definitively state that it was from the phase transition we 
considered. However, there have been extensive investigations of existing weak-scale models, 
primarily with the aim of studying electroweak baryogenesis, and no model has yet been 
found with such a strong transition (i.e. with a ^> 1). For example, the phase transition we 
have considered can be much stronger than the usual EW phase transition, which is usually 
constrained by the experimental bound on the Higgs mass. In our case, there are no such 
bounds. The radion sector is much less constrained and this offers the possibility of a strong 
phase transition. Although it would not be conclusive, a large LISA signal of the sort we 
have discussed is likely to be associated with an RSI cosmo logical phase transition. 

This gravitational wave signature has the advantage of being common to all RSI models. 
This is to be contrasted with collider signatures which depend on the details of the model 
such as the localization of the Standard Model fermions and gauge fields in the bulk of AdS§. 
On the other hand, the signal strongly depends on the radion sector which to some extent 
can be probed by collider experiments. We are finding that the phase transition can proceed 
only in the regime of large back reaction. According to Ref. [24], the radion couplings are 
expected to be large in this case. This could indeed imply the possibility of testing the 
strength of the transition at colliders. 
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We emphasize that the uncertainties in the calculation are primarily due to the proximity 
to the perturbativity limits of the calculation. However, if the phase transition completes, it 
is clear that the signal is likely to be large. It will be interesting to see how far these results 
can be extended into the nonperturbative regime. It will also be interesting to investigate 
other possibilities for weak scale physics and gravitational waves. 
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